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ABSTRACT 

Almost all of the extragalactic X-ray background (XRB) at 0.25 keV can be accounted for by 
radio-quiet quasars, allowing us to derive an upper limit of 4 keV cm~^ s~^ sr~^ keV~^ for 
the remaining background at 0.25 keV. However, the XRB from the gas halos of groups of 
galaxies, with gas removal due to cooling accounted for, exceeds this upper limit by an order 
of magnitude if non-gravitational heating is not included. We calculate this using simulations 
of halo merger trees and realistic gas density profiles, which we require to reproduce the 
observed gas fractions and abundances of X-ray clusters. In addition, we find that the entire 
mass range of groups, from '^ 5 x 10^^ to ^ 10^'^Mq, contributes to the 0.25 keV background 
in this case. 

In a further study, we reduce the luminosities of groups by maximally heating their gas 
halos while maintaining the same gas fractions. This only reduces the XRB by a factor of 2 
or less. We thus argue that most of the gas associated with groups must be outside their virial 
radii. This conclusion is supported by X-ray studies of individual groups. 

The properties of both groups and X-ray clusters can be naturally explained by a model in 
which the gas is given excess specific energies of ^^ 1 keV/particle by non-gravitational heat- 
ing. With this excess energy, the gas is gravitationally unbound from groups, but recoUapses 
with the formation of a cluster of temperature > 1 keV. This is similar to a model proposed 
by Pen, but is contrary to the evolution of baryons described by Cen & Ostriker. 

The excess energy is most likely injected by galaxies in the smaller 'branches' of the halo 
merger tree (< 5 x IO^^Mq), by active galactic nuclei and possibly supernovae. The heating 
process may therefore play an important part in the evolution of galaxies. 

In addition to the soft XRB spectrum, we simulate source counts in two bands: 0.1-0.4 
keV and 0.5-2 keV, for comparison with present and future data. 

Key words: galaxies: clusters: general - galaxies: haloes - intergalactic medium - cooling 
flows - X-rays: galaxies 



1 INTRODUCTION 

Compared to X-ray clusters, relatively little is known about the hot 
gas halos of galaxy groups (henceforth simply 'groups'). With the 
Chandra and XMM satellite missions we can expect much to be 
revealed about them. However, the extragalactic soft X-ray back- 
ground (XRB) below ~ 1 keV already provides a useful probe 
of their mean properties. For example, we shall show that (in the 
absence of non-gravitational heating) the 0.25 keV background 
probes almost the entire mass range of groups (~ 5 x 10^^- 
10^'' M0, corresponding to temperatures of T ~ 0.1-1 keV). 

Groups are cosmologically important, for a majority of the 
bary ons in the universe may be associate d with them at the present 
day (Fukugita, Hogan, & Peebles 1998). In addition, groups con- 



nect galaxies to clusters in a hierarchical merger tree. This suggests 
that a large fraction of the heavy elements and 'excess energy' 
(due to non-gravitational heating) injected into the inter-galactic 
medium (IGM) would have to pass through groups before ending 
up in X-ray clusters. However, it has been shown that strong heat- 
i ng of the IGM is required to satisfy constraints on the soft XRB 
(Pen 1999). The heating required by Pen, 5T ~ 1 keV, would im- 



ply that most of the gas would not be gravitationally bound to the 
potential wells of groups, contrary to the scenario described above. 
It is also interesting that this amount of heating is similar to that 
required to e xplain the properties of X-ray clusters (W u, Fabian, & 
Nuls en 1998';^oewenstein 1999|; |Pen 1999|; W u, Fabian, & Nulsen 
1999). Together they therefore suggest a consistent picture of heat- 
ing in the IGM. Furthermore, the mechanism for heating the IGM 
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is likely to have significant implications for the evolution of galax- 
ies (Wu et al., 1999; hereafter WFN99), whether the heat source be 
supemovae or active galactic nuclei (AGN). 

In this paper, we examine in more detail the contraints on 
groups, using simulations of halo merger trees and drawing on el- 
ements of our semi-analytic model of galaxy formation (WFN99). 
One important difference from Pen's model is our inclusion of cool- 
ing, in particular its role in removing gas from halos. Our methods 
for calculating the XRB also differ significantly. Broadly speaking, 
the main questions we aim to address are whether groups are signif- 
icantly affected by non-gravitational heating, and if so, what range 
of groups are affected. 

Although a large fraction of the soft XRB below about 0.8 
keV originates in our own Galaxy, the extragalactic component can 
be measured with sh adowing experiments. As reviewed by W ar- 
wick & Roberts (1998), measurements of the extragalactic XRB in 
the 0.1-0.4 keV band (hereafter 0.25 keV) appear to be converg- 
ing on a value in the range 20-35 keV cm~^ s~^ sr~^ keV~^. 
However, the bulk of the extragalactic XRB is due to AGN. A sig- 
nificant fraction of the XRB at 1 keV has been resolved into AGN. 
By r nodelling the OSO X-ray luminosity function and its evolu- 
tion, Boyle et al. (1994) e stimated the contribu tion of QSOs to the 
1-2 keV XRB (see also [Schmidt et al. 1998| ). Their results cor- 
res pond to 35-55 per cen t of th e 1-2 keV backg round measured 
by 3endreau et al. (1995 ). Now, Laor et al. (1997 ) find that radio- 
quiet quasars have a mean spectral slope of a = —1.72 ± 0.09 
over the range 0.2-2 keV, though they obtain a significantly flatter 
slope for (less common) radio-loud quasars. Therefore, if we as- 
sume that radio-quiet quasars account for 90 per cent of the quasar 
contribution to the 1 keV background, then they contribute at least 
0.9 X 0.35 or about 30 per cent of the 1 keV background . Using 
the 1 keV background measured by Gendreau et al. (1995), which 



is 9.6 keV cm^ s^ sr~ keV^ , we can thus obtain a lower limit 
to the contribution of radio-quiet quasars below 1 keV. At 0.25 keV 
the lower limit is 31 keV cm^^ s^^ sr~^ keV~^. From the range 
quoted by Warwick & Roberts, it follows that gas halos almost cer- 
tainly contribute less than 4 keV cm~^ s~^ sr~^ keV^^at 0.25 keV, 
as most of the 0.25 keV background is already accounted for by 
QSOs. We shall adopt this value as an upper limit to the contri- 
bution from gas halos. (We note that Boyle et al., Laor et al. and 
the detections quoted by Warwick & Roberts all obtained their data 
from the same instrument, namely the ROSAT Position Sensitive 
Proportional Counter.) The main assumption we have made lies in 
the application of the mean properties measured by Laor et al. to all 
quasars. On the other hand, we have used the lowest expected con- 
tribution of QSOs at 1 keV, along with the highest expected value 
for the 0.25 keV background. We do not attempt here to assign a 
confidence level to this limit since the main uncertainties are sys- 
tematic and not statistical. 

The models in this paper are intended to be realistic, but sim- 
ple enough for comparisons to be easily made. We consider two 
low-density universes with $!,„ = 0.3: one with no cosmologi- 
cal constant and the other with JIa = 0.7. A Hubble constant of 
Ho = lOOh km s~^ Mpc^^ where h = 0.5 is assumed (the in- 



sensitivity of our results to h is discussed in section 3.1). We use 



CDM-type fluctuation power spectra, normalized to match the ob- 
served abundance of X-ray clusters. In this way, the abundance of 
groups are in effect extrapolated from the X-ray cluster abundance. 
Likewise, in our fiducial models we assume that groups have a gas 
fraction of 0.06ft~'^ " when they virialize, as this is the mean value 
measured for clusters (Evrard 1997; Ettori & Fabian 1999). This 



Gas halos are assumed to be isothermal and hydrostatically sup- 
ported in Navarro, Frenk & White (1997; NFW) potential wells. 
The resulting gas profiles are found to closely model X-r ay clus- 
ters (iMakino, Sasaki, & Suto 19971; lEttori & Fabian 19991). In the 
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absence of excess energy, the model gas profiles of groups and 
clusters are self-similar to a good approximation. This is a reason- 
able assumption to make (see hydrodynamic simulations such as 
Navarro, Frenk & White 1995), especially as we are looking for 
large differences between model results and observations. We syn- 
thesi ze spectra using the MEKAL spectral synthesis code (K aastra 
1992). 

This paper is organised as follows: in sections we describe our 
model, followed by how we calculate the XRB and source counts. 
We simulate the XRB from 0.05-2 keV, and calculate source counts 
in two bands: 0.1-0.4 keV and 0.5-2 keV. In most of sections Q and 
we assume that non-gravitational heating is absent. However, we 
also describe simulations where halos are required to follow ob- 
served luminosity-temperature relations extrapolated to lower tem- 
peratures. The fiducial results are discussed in section H. In sec- 
tion H we investigate the effects of including non-gravitational heat- 
ing, and of allowing the gas fraction of halos to be determined by 
inheritance. We discuss some implications of our results in sec- 
tion N and summarize our conclusions in section o. 



2 SIMULATION 

The main components in the calculation of our fiducial results are 
the halo merger tree, the gas density profiles of virialized halos, and 
the spectral synthesis model. 

The halo merger trees are simulated using the C ole & Kaiser 
(1988) block model, as in our earlier models (e.g. WFN99). The 
smallest regions simulated have a mass of 1.5 x 10^" Mq. We as- 
sume a collapse hierarchy of 20 levels, so that the mass of the 
largest block is 2^'^ x 1.5 x lO"' = 7.9 x W^^Mq. This is the 
total mass of the region simulated in one 'realisation' of the merger 
tree. 

The mean density of collapsed halos are specified to conform 
exactly with the spherical collapse model. In open cosmologies 
without a cosmological constant, the mean density of a virialized 
halo is assumed to be IStt^ ~ 178 times the background density of 
an Einstein-de-Sitter universe of the same age. In flat cosmologies 
with a cosmological constant, we use the analytic approximation 
given by Kitayama & Suto (1996; equation A6), where the mean 
density of a halo is given by 



Pvir = 187r^(l + 0.4093(n„^ - 1)° 



')pb 



(1) 



is a reasonable choice since groups are the progenitors of clusters. 



Here, the density parameter, flm, and the background density, Pt, 
are evaluated at the time of collapse. At very early times, results 
from the two prescriptions tend to the same value, as the vac- 
uum density is then small relative to the density of collapsing ha- 
los. However, at late times the simple prescription for open cos- 
mologies becomes a poor approximation in A-cosmologies. For 
example, for a halo that virializes today in a cosmology given 
by {Q,rn,^A,h) = (0.3,0.7,0.5), the first prescription underes- 
timates the mean density by about 20 per cent. Since the X-ray 
luminosity of a gas halo scales roughly as the density squared, this 
difference can be significant. From the mass of the collapsed halo, 
the radius of the halo, rvir, thus follows. 

The second component in the calculation is the assumed gas 
density profile of halos. To simplify the synthesis of spectra, we 
consider only isothermal gas halos. We assume the gas to be in 
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hydrostatic equilibrium in Mavarro, Frenk, & Wliite (1997; NFW) 



potential wells. In other words, the total density of the halo is de- 
scribed by the NFW profile: 



p{r) = 



ScPc 



(r/r.)(l+r-/r,)2' 



(2) 



where pcrit = SH^ /SttG and H is the Hubble parameter at the 
time of collapse. The parameter Sc is calculated as described in the 
Appendix of NFW. The scale radius r, is then uniquely determined 
by the mean density calculated above. The concentration parameter, 
c, is defined to be rs/rvir. (This differs slightly from NFW, as we 
use a more detailed derivation for the mean density. In particular, 
their equation (2) that relates Sc and c is slightly modified in this 
model.) From the NFW profile, the gravitational potential is given 
by 



Uc = a - 



ln(l + a;) 



1 + c 



(3) 



where x — r/rs and a = A-rGpsTs ■ It then follows that a gas halo 
of temperature T in hydrostatic equilibrium takes the form 



Pg^ll+x)"/^ 



(4) 



where pg is the gas density, and r) — ^vana/ikT) (Wu et al., 
1998; WFN98). Here, the mean mass per particle of the gas is de- 
noted by /iTTiH, and k is the Boltzmann constant. This gas pro- 
file closely approximates the co nventional /3-model if /3 = '7/15 
( Making, Sasaki, & Suto 1997 ), and models most X-ray clusters 
very well. For large clusters the m ean value of r) is observed to be 
about 10.5 dettori & Fabian 1999| ). 

We calculate rj by first specifying the total energy of the gas 
halo, which then uniquely determines t). This is described in more 
detail elsewhere (WFN98; Model B in WFN99). To summarize, in 
the absence of excess energy, the total specific energy of the gas 
halo (thermal plus gravitational) is required to be proportional to 
the specific gravitational energy of the entire halo; we then cali- 
brate this energy relation by matching to the largest clusters. Al- 
though halos described by the NFW profile are not exactly self- 
similar (as the concentration c varies), this relation for the gas halos 
expresses one form of self-similarity. For the purposes of this pa- 
per, the main point to note is that this results in r] close to 10.5 for 
all gas halos. In addition, the resulting gas temperatures are well- 
approximated by Tsis, the temperature that the gas would have if 
both gas and dark matter had power-law density profiles: p on r^'^; 
i.e. kTsis/{pmu) — GAftot/(2rvir), where Mtot is the total mass 
of the halo. In general, T scatters between Tsis and l.OSTsis for 
clusters, with the upper end of the range increasing to l.lTgig as 
we go down to halo of ~ 10^^ Mq. 

We assume that isolated galaxies have halo masses of up to 
~ 10^^ Mq. Since we are primarily concerned with the properties 
of halos, in this paper the term 'group' simply refers to a halo of 
mass greater than a few W^^Mq (but less than ~ 10^'' Mq). 

In our fiducial models, we assume that all halos have a gas 
mass fraction of 0.17 wh en they viriali ze. This is the mean cluster 
gas fraction measured by Evrard (1997) and Ettori & Fabian (1999) 
(using h — 0.5), within a radius rsoo defined to be such that the 
mean density within rsoo is 500pcrit . We note that the gas fraction 
of both observed and simulated clusters tend to increase with the 
radius within which they are measured. Since rsoo < '"vir, it is 
therefore likely that the true gas fractions of clusters within r^ir are 
higher than this. This strengthens our argument below as it would 
increase our fiducial estimates of the XRB. 



2.1 Calculating the X-ray background 

We now describe in detail how we calculate a mean spectrum for 
the X-r ay backgroun d. Spectra are simulated using the MEKAL 
model (Kaastra 1992), to which we supply two parameters: the gas 



temperatur e, T, and the metallicity, Z, in units of the solar abun- 
dance, Zq (Anders & Grevesse 1989). From the spectrum given by 



MEKAL, we calculate the cooling function, A(r, Z). In this way, 
the model self-consistently estimates the cooling time of the gas. 



2.1.1 Cold and hot collapses 

We recall that in the collapse of less massive halos, the cooling time 
of the gas, tcooh can be shorter than the free-fall time to the centre 
of the halo, iff. In this case the gas cools fast enough that it is not 
hydrostatically supported, and in spite of any shock heating, the gas 
temperature remains well below the virial temperature most of the 
time. We refer to this as a cold collapse (WFN99). In our model, 
gas is labelled as cold when r = tcooi/iff < to, where tq ^ 1 
is a free parameter. Otherwise, the gas is labelled as hot. This is 
important when calculating the X-ray background, as we expect 
cold gas to radiate little, if at all, in the X-ray band. Therefore we 
only integrate contributions from hot gas and use a fiducial value 
of To = 1. 

For isothermal gas profiles, r is almost always a monotoni- 
cally increasing function of radius (see Appendix A of WFN99 for 
a more detailed discussion). This allows us to define a radius, rcp, 
such that gas outside rcp is labelled as hot and gas inside as cold. 
As halo mass increases, rcp moves from outside the virial radius 
to the centre of the halo. This transition is fairly abrupt and occurs 
over about one decade in mass, in halos of ~ lO^^M©. For ha- 
los in the transition region, rcp is found by solving the equation 
''"('"cp) = To numerically. 

(In our model, the transition to the hot-gas regime is made 
more abrupt by supernova feedback from star formation. We as- 
sume that cold gas rapidly forms stars, which lead quickly to type 
II supernovae. If a sufficient fraction of a collapse is cold, then the 
energy from supernova feedback is able to eject the res t of the at- 
mosphere, including the hot gas (Nulsen & Fabian 1997). By using 



a much lower value of to = 0.1, we show that these complications 
at the transition region have a very small effect on the predicted 
XRB and source counts, and do not affect our conclusions in any 
way.) 

When a hydrostatically-supported hot gas halo occurs, any gas 
that cools is assumed to form low-mass stars or 'baryonic dark mat- 
ter', in analogy to cooling flows in X-ray clusters. A possible mech- 
anism for low-mass star form ation in cooling flows is described by 
Mathews & Brighenti (1999), for the case of elliptical galaxies. It 



remains possible (if not likely) that some normal star formation and 
feedback occurs in cooling flows. However this does not affect our 
main conclusions (the effects of strong heating are investigated in 
section ft). 



2.1.2 Spectral synthesis 

We calculate model XRB spectra from E = 0.05 keV to 2 keV. 
We first divide this range into equal logarithmic bins of width 
A logj^Q E — 0.01. Suppose photons in the universe belonging 
to the energy bin {E, E + AE) at the present day have a num- 
ber density of An, then the corresponding energy flux per stera- 
dian (i.e. the intensity) is given by EAnc/{4n), where c is the 
speed of light (in practice, we use the geometric mean of E and 
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radius 



Figure 1. Schematic diagram of liow the total energy radiated by a halo is 
computed. The solid curve encloses the gas emitting in X-rays at any given 
time: between rcp and TcooI it gives the cooling time of the gas, but gas 
outside Tcool does not have time to cool before the next collapse, which 
occurs at tn- Gas inside rcp is labelled as cold and therefore does not 
contribute to the XRB. The correspondence between rest-frame bins and 
collecting bins shifts by one position at the times ti, . . . , tn-i- Thus the 
energy radiated during each time interval needs to be calculated. For the 
gas between tqy ^nd rcool . our method involves first dividing the halo into 
shells, as shown by the dashed lines (see text). 



E + Ai5 in place of E in this expression). Dividing by Ai? thus 
gives the intensity per unit energy, which we express in units of 
keV cm^^ s^^ sr^^ keV^^. Below, we refer to these bins as 'col- 
lecting bins'. 

Since each realisation of a merger tree simulates a region of 
constant comoving volume, the photon density An is simply given 
by the total number of photons (of the correct energy) emitted in the 
simulation divided by the present-day volume of the simulation. 

We now consider a gas halo at some redshift z. Given T and 
Z, its rest-frame spectrum is calculated with MEKAL, also in equal 
logarithmic bins of width A logj^Q E = 0.01. The spectrum is inte- 
grated to give the cooling function A(r, Z) — defined such that the 
bolometric luminosity per unit volume is given by n^nYiA.{T, Z), 
where n^ and nn are the electron and hydrogen number densities 
respectively. We then calculate the number of photons emitted in 
each 'rest-frame bin' per unit total energy that is radiated. These 
ratios are used throughout the remaining calculation. 

For gas radiating at a redshift of z, each collecting bin is 
blueshifted accordingly (from E lo {1 + z)E) to collect pho- 
tons of the correct energy. Photons from each rest-frame bin are 
then dumped into the nearest collecting bin. Since the amount of 
blueshift is a continuous function of z, the correspondence be- 
tween rest-frame bins and collecting bins shifts by one bin each 
time (1 + z) decreases by 0.01 dex. This occurs ~ 10 times over 
the life of a halo (which is the period from one collapse to the next 
in the block model — we interprete collapses as major mergers). 

To summarize, for each halo we calculate the total energy ra- 
diated during each time interval between bin shifts. From this the 
number of photons deposited in each collecting bin is found. The 
total number of photons collected in the simulation then gives the 
XRB spectrum as explained above. 

In what follows we describe how we calculate the energy radi- 
ated by a halo, taking into account the effects of cooling. The most 
general case is illustrated in Fig. hi which shows the amount of gas 



contributing to the XRB as a function of time. Also shown are the 
times tiii = 1,2, ... , n~ 1) when bin shifts occur, measured from 
when the halo virializes. t„ gives the time of the next collapse. 

Gas inside rcp is regarded as cold and therefore is not in- 
cluded in the calculation. Between rcp and the cooling radius, 
rcooi, the cooling time is shorter than the time to the next collapse: 
gas in this region is assumed to cool out once its cooling time has 
elapsed. The cooling time is estimated by 



tc 



3 pgkT/fimji 
2n~mJMT\Z)' 



(5) 



As explained in WFN99, the gas profile that we use to estimate 
such quantities must be regarded as notional, as it cannot describe 
the gas halo at all times. In particular, the gas halo redistributes 
itself as gas in the centre cools out. Nevertheless, the gas profile 
allows us to estimate the behaviour of different subsets of gas in 
the halo. 

Between rcooi and r^ii, gas does not have time to cool before 
the next collapse. Here, the energy radiated during each time in- 
terval is easily calculated from its bolometric luminosity, which is 
given by J "" nenHA(r, Z)dV. The number of photons dumped 
into each collecting bin thus follows. 

This is less straightforward to apply to the gas inside rcooi, be- 
cause the amount of gas changes continuously with time. However, 
the gas inside rcooi cools completely, so that the energy radiated 
is simply given by {3/2)NkT, where A'^ = J pg/{iimii)dV is 
the number of particles that cool. To estimate the amount of energy 
radiated in each time interval, we first divide the halo into shells 
separated by the radii Rcooi{ti). The function -Rcooi (i) gives the 
radius where cooling time is equal to i. For the shell with an outer 
radius of -Rcooi (ij)' ws assume that {tj — ij_i)/fi of the energy 
it radiates (as given above) is emitted in the j-th time interval. In 
this way, the total energy radiated is correctly accounted for, but 
the 'allocation' of energies between time intervals is less precise. 
However, the latter is equivalent to determining how the photons 
are binned, and the binning of photons is already uncertain by ±1 
bin, because the bin shifts are uncontinuous. It follows that the de- 
rived spectrum is 'smeared' by ~ 1 bin width. 

Note that the release of (3/2) fcT per particle that cools must 
be regarded as a lower limit: in reality, the weight of the overly- 
ing gas is likely to at least maintain the pressure of the gas as it 
cools, so that gravitational work raises the total energy radiated to 
(5/2) fcT per particle or more in most cases. In addition, the den- 
sity of the gas increases as it cools and moves to smaller radii, so 
that the luminosity of a given parcel of gas increases. As a result, 
(3/2) fcr leads to a reasonable estimate of the cooling time. For 
the gas outside rcooi, the density also increases as the gas moves 
inwards. Therefore, our estimate of the energy released in this case 
should also be a lower limit. 

A caveat of our model is that any intrinsic absorption that 
might occur in the cooling region (in analogy to cluster cooling 
flows) is ignored. Since this is confined to the cooling region at any 
given time, while most of the simulated XRB arises from outside 
that region, its effect is unlikely to be important. 



2.2 Source counts 

The log -/V-log S function for X-ray halos is calculated in two 
bands: the 0.1-0.4 keV band, for comparison with future results 
from Chandra and XMM, and the 0.5-2 keV band, for comparison 
with results from the Wide-Angle ROSAT Pointed X-ray Survey 
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( WARPS; Scharf et a l. 1997) and the slightly deeper counts made 
by|Rosatietal. (1995|). 



Since the simulated merger trees provide no spatial informa- 
tion, it is assumed that the halos are distributed randomly in space. 
The probability of a given source in the simulation being observ- 
able in principle (with an infinitely sensitive telescope) depends on 
how long it 'exists'. Suppose a source exists for a short period cor- 
responding to dz. The comoving volume observable by us on the 
entire sky in this redshift interval is given by 



dVe 



[D^{z)/{l + z)f 



Ho y/nm{i + zyi + (1 - n„, - nA)(i + zy + Qa 



inldz 



where Dl (z) is the luminosity distance (given below). Suppose we 
create an infinite universe by tiling together copies of the same sim- 
ulation. Then the mean number of copies of this particular source 
that we expect to see is given by dVc divided by the volume of sim- 
ulation. This is therefore the contribution of this source to the mean 
source count on the whole sky. 

We integrate source counts for fluxes of S = 10~^^- 
10"^"^ erg s~^ cm~^, which we divide into logarithmic bins of 
width Alogjo S — 0.1. The flux from a gas halo is given by 
S = I/x/(47rZ)|), where Lx is the luminosity in the relevant, 
blueshifted band. For JIa = 0, Dl is given by (e.g. Peacock 1999): 



Dl 



(1 + 2 



but for Q„ 



^L = (i + ^)Tr 



2c D,„,z + {n^ ~ 2)(V1 + QmZ - 1) 
'Hq n^l + z) 

Qa = 1, Dl is integrated numerically: 

dz' 



Ho Jo ^n^{i + zy^ + nA 



(V) 



(8) 



Equation |6| is integrated over the lifetime of a given halo in the 
simulation, to give the mean number of copies of that halo that we 
expect to see on the sky. To do this we associate Az with the same 
time intervals described in the previous section. This gives a good 
approximation for the cosmologies and redshifts that we are inter- 
ested in. Luminosities and fluxes are calculated at the beginning 
of each time interval. Each interval then contributes to the num- 
ber count in the relevant flux bin as described above. (Using the 
average luminosity in each time interval only changes the number 
counts by ^ 5 per cent.) 

We tested our code by summing all the discrete sources to give 
an alternative calculation of the XRB at a given energy. To do this 
we replaced the above energy bands with one very narrow band at 
this energy, and integrated the total flux in this band. The resulting 
intensity per keV agreed well with the spectrum derived using our 
usual method. 



2.3 Empirical calculations based on Lx T relations 

It is well-known that self-similar gas halos (such as those described 
at the beginning of this section) do not match the observed proper- 
ties of X-ray clusters. Substantial heating is required to lower the 
luminosities of the smaller clusters, in order to match the observed 
relation of Lx oc T^ (WFN98,99; Pen 99). How this extends to 
groups is less clear. A clue is offered by the Lx — T distribution 



than in clusters. Other X-ray studies of groups ( Mulchaey et al. 
1 996; [Helsdon & Ponman 1999[ ) confirm this view. 

Since we only require the temperature and bolometric lumi- 
nosity of a gas halo in order to estimate its spectrum, we have per- 
formed separate simulations based on (extrapolations of) observed 
Lx — T relations. A basic difficulty with the Lx — T relation for 
groups is that a large fraction of groups do not have detectable X- 
ray emission, so that an X-ray selected sample is likely to give a bi- 
ased estimate of the average luminosity. In fact, optically-selected 
samples also face difficulties, as they may include chance projec- 
,(6)tions of galaxies. On top of this, the X-ray temperature of a group 
can differ significantly from the temp erature associated w ith the 
velocity dispersion (Tofthegalaxi£s(Ponmaneta^ 
don & Ponman 1999). If we suppose that a reflects the underlying 
dark matter distribution, then the gas temperature in the self-similar 
case is approximately given by kT/{fimu) ~ cr^. In reality, the 
T — a distribution for groups exhibits a large scatter about the best- 
fit power law, which has a different slope from the self-similar re- 
lation. For the HCGs, the best-fit power law yields temperatures 
(as a function of a) th an are up to a facto r of 2 higher than the 
self-similar prediction ( Ponman et al. 1996| ). 

Faced with these uncertainties, we have simply applied the ob- 
served Lx — T relations and assumed that T is the same as the tem- 
perature given in the self-similar model. In addition, we assume that 
the luminosity of a halo remains constant during its lifetime. This 
is reasonable because the time-averaged luminosity of a halo in 
the simulation should correspond to the observed luminosities av- 
eraged over many groups. In this way, we are able to calculate XRB 
spectra and source counts as described above. Although crude, the 
results give an insight into what the true composition of the soft 
XRB might look like. 

The best-fit Lx — r relation for the HCGs is given by (P onman 
et al. 1996): 



of the Hickson Compact Groups or HCGs (Ponman et al. 1996), 



which is significantly steeper than that for clusters. However, this 
only extends down to about 0.5 keV, and the scatter in the distri- 
bution appears to be intrinsically very large. (In section W we pro- 
pose a reason for the large scatter.) The steeper Lx — T relation 
for groups suggests that the effects of heating are even more severe 



log Lb 



(43.17 ± 0.26) + (8.2 ± 2.7) log T, 



(9) 



where Lboi is the bolometric luminosity in erg s~^ and T is in 
keV. As noted above, we also apply this relation below its ob- 
served range. Despite the uncertainty in the slope, it is clearly much 
steeper t han the relation for clusters. For c lusters, we use the best-fit 
given by White, Jones, & Forman (1997): 



log Lb 



(42.7 ± 0.1) + (2.98 ± 0.11) log T. 



(10) 



This spans a range of T ~ 2-10 keV. The two relations intersect at 
r = 0.8 keV. Thus we are able to apply the cluster relation above 
0.8 keV, and the group relation for lower temperatures. We refer to 
this combination as 'PW'. We also perform simulations with only 
the cluster relation, extrapolated to all temperatures. We refer to 
these as 'W' simulations. 

In the simulations, we do not include scatter in the above re- 
lations. This is permissible for estimating the XRB, which depends 
on the mean properties of halos. However scatter can significantly 
affect the resulting source counts. For our purposes, this is more 
easily adjusted for after the simulation, and we discuss it with the 
results. 



3 FIDUCIAL RESULTS 

3.1 Calibration 

Perhaps the two most important 'parameters' in this study are the 
gas fraction and the mass function of the relevant halos. 
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Table 1. List of simulations and legend for Figs. wA and Figs. H-llll Each 
line in the table corresponds to two simulations, performed in the OCDM 
and ACDM cosmologies. Only the parameters that change between simu- 
lations hav e been listed. The 'W and 'PW simulations are described in 
section |2.3[ 

line-style metallicity (Zq) tq F Lx — ^ fixed? 



sohd 


0.3 


1 


0.106 


no 


dotted 


0.03 


1 


0.106 


no 


dashed 


0.3 


1 


0.25 


no 


dot dash 


0.3 


0.1 


0.106 


no 


dash dot dot dot 


0.3 


1 


0.106 


yes,W 


long dashes 


0.3 


1 


0.106 


yes, PW 



For the fiducial models discussed in this section, we set the 
initial gas fraction of all halos to be 0.17. Although this may be 
unrealistic in detail, it provides a simple point of reference. As for 
the mass function of groups, the primary constraint, albeit indirect, 
is the mass function of X-ray clusters at 2 = 0. We require all 
of our simulations to reproduce the temperature function of X-ray 
clusters, fo r which we use the best-fit power law measured by E dge 
et al. (1990). (The temperature function is a reasonably direct way 
of constraining the mass function of clusters, whereas the luminos- 
ity function is sensitive to other properties, such as the gas fraction 
and the structure of gas halos.) It follows that the abundance of less 
massive halos are in effect extrapolated from the observed cluster 
abundance, by assuming CDM power spectra. To investigate the in- 
herent uncertainty in this, we consider two different values for the 
CDM shape parameter V. 

As in WFN99, w e assume a baryon density parameter o f 
0.08 (|Burles & Tytler 199^ : 



Buries etal. 1999) 



fib = 0.02h 

Following jugiyama (1995| ), the CDM shape parameter is then 
given by r = flmhexp{-flb{l + V2h/Qm)) = 0.106. Using this 
value ofr in an OCDM universe with (n^.f^ A, ft) = (0.3,0,0.5), 
we are able to match the temperature function of X-ray clusters us- 
ing o-g = 0.75 (as shown in WFN99). 

For a ACDM cosmology with (fi^, D.A,h) = (0.3, 0.7, 0.5) 
and the same values for Qi, and F, we are able to reproduce the 
cluster temperature function using ag = 0.8. Both of the above 
results for ag agree remarkably well with Eke et al. (1998, and 
private communication for ACDM). 

For each cosmol ogy, we also consider a hi gher value of F = 
0.25, as measured by Peacock & Dodds (1994) from galaxy sur- 
veys. In both cosmologies, we find that the same values of ag are 
able to reproduce the cluster temperature function, in agreement 
with Eke etal. (1998). 

A list of simulations is given in Table 111 which we repeat for 
each cosmology. The table gives the values of parameters that vary 
between simulations. We assume the same metallicity for all gas 
halos in any given simulation. For each simulation that assumes 
no non-gravitational heating ('no' in the last column), we use 100 
realisations of the merger tree. These simulations produce similar 
results to each other. Thus, to test that their results have converged, 
we increased the number of realisations in one of these simulations 
by a factor of 4. This made no noticable difference to the results. 

In the simulations based on extrapolations of observed Lx — 
T relations ('yes' in the last column), the XRB is dominated by 
more massive halos, so that 400 realisations are required to obtain 
convergence. 

Finally, we note that our results are not sensitive to increases 
in h. For our fiducial simulations that assume self-similarity, this 
can be understood by noting that the model is normalized to repro- 
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Figure 2. Simulated XRB spectra in the OCDM cosmology. The legend 
for the different lines is given in Table R The upper limit for the gas halo 
contribution to the observed 0.25 keV background is shown by the arrow. 



duce the contributio n of the largest obser ved clusters (which are 
close to self-similar, A.llen & Fabian 1998) to the XRB. Since that 



contribution is a directly measured quantity, it has no dependence 
on h. The rest of the simulated XRB is then a non-trivial extrapo- 
lation of the cluster contribution. A similar argument applies to all 
our simulations. 

To test the sensitivity, we also simulated the first model in Ta- 
ble \n using h — 0.7, without changing F and using ACDM. The 
gas fraction of clusters scales as h^"^'^, but this is compensated 
by the number density of halos — as expressed by the temperature 
function — which scales approximately as h'^. The resulting XRB 
spectrum is about 10 per cent higher and the source counts (in the 
given flux range) are virtually unchanged. 



3.2 The XRB spectrum from hot gas halos 

The XRB spectra from the simulations listed in Table hi are dis- 
played in Figs. H and ||, for OCDM and ACDM respectively. 
It is clear from the two figures that the two cosmologies pro- 
duce similar results. (The OCDM spectra are only slightly higher 
than the ACDM spectra, by up to a factor of 1/3 depending 
on the energy.) The 0.25 keV backgrounds predicted by the 
simulations that assume no heating are all much higher than 
our upper limit of 4 keV cm^^ s^^ sr^^ keV^^ (section 1) — 
several by over an order of magnitude. These simulations also 
overpredict the I keV background, which should be about 1 
keV cm~^ s~^ sr~^ keV~^; the extragalac tic background at 1 keV 
is about 10 keV cm"^ s"^ sr"^ keV"^ (pendreau et al. 19951), 



but only about 10 per cent of this is likely to be from groups and 
clusters (iMcHardy et al. 19981). 



To illustrate the importance of cooling, we modified one of 
these simulations as follows to switch off cooling. For all halos that 
contribute to the XRB in the fiducial simulation, a) we label all of 
their gas as hot (thus rcp = 0), and b) we suppose that no gas is 
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Energy (keV) 
Figure 3. As Fig. Q but in the ACDM cosmology. 



removed as a result of cooling. This increases the 0.25 keV back- 
ground by a factor of 5, and increases the 1 keV background by 
almost 3-fold. These large increases show that if cooling is not in- 
cluded, then gas will often end up radiating many times its thermal 
energy. (In this case, our results increase to a level comparable to 
the 0.25 keV background calculated by Pen (1999), using clumping 
factors obtained from hydrodynamic simulations.) 

We shall now discuss the differences between the above spec- 
tra. We regarded the parameters used to obtain the solid spectra as 
the 'default' parameters (Table ll|) and varied each of the parame- 
ters in turn to observe the consequences. Reducing the metallicity 
from 0.3 to 0.03 Zq (dotted spectra) reduced the amount of line- 
emission, so that the resulting XRB spectrum is much smoother. 
By comparison, a large 'bump' in the solid spectrum at around 0.7 
keV can be seen. It corresponds to the (redshifted) iron L com- 
plex. However, the changes at the other energies are more modest; 
the spectrum is almost independent of metallicity at 0.25 keV. It 
is interesting that the XRB spectrum increases at certain energies 
when we reduce the metallicity. This occurs in parts of the spec- 
trum where line-emission is weak. It is possible because the cool- 
ing time of gas increases when we reduce the metallicity. Increas- 
ing the CDM shape parameter to F = 0.25 (dashed spectra) raises 
the power at sub-cluster scales, resulting in more galaxy groups. In 
both cosmologies, this raises the XRB spectrum by roughly 50 per 
cent. 

To better understand the implications of these results, we need 
to know which types of halo are contributing to the soft XRB. We 
therefore computed the fractional contribution of halos belonging 
to each mass in the block model. This was done at two energies: 
0.25 keV and 1 keV, which correspond to the two energy bands 
used in our source counts. In Figs. H and dI we show the fractional 
contributions to the solid spectrum in the ACDM case (Fig. |3|), at 
0.25 keV and 1 keV respectively. The corresponding plots for the 
OCDM case are almost the same, and the results for the dotted and 
the dashed spectra are very similar. In these plots, the solid lines 
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Figure 4. Plot of halo contributions to the XRB at 0.25 keV, for the case 
of the solid spectrum in the ACDM cosmology (Fig. pi). The solid lines 
correspond to the left axis, and the dotted line to the right axis. The squares 
give the fractional contribution from all halos of a given mass, and the trian- 
gles give the mean contribution from those halos. The squares in effect give 
the contribution per unit logarithmic interval in mass. The crosses show the 
number of hot gas halos obtained in the simulation; thus the squares are the 
product of the triangles and crosses. 
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Figure 5. As Fig. W but for the XRB at I keV. 
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Figure 6. As Fig. W but for the dot-dashed spectrum at 0.25 keV in the 
ACDM cosmology. Lowering to to 0.1 results in hot gas halos at lower 
masses. In addition, the low-mass drop-off in the squares is now due to 
the sharp decline in contribution per halo (shown by the triangles), which 
should be contrasted with Fig. W. 



correspond to the axis on the left and the dotted line to the axis 
on the right. The main result is given by the squares, which show 
the fraction of the XRB due to all halos of a given mass (since the 
masses increase by factors of 2, this is equivalent to the contribution 
per unit logarithmic interval in mass). 

At 0.25 keV (Fig. W), the squares form a relatively flat plateau 
from ~ 5 X IO^^Mq to ~ 10^^ Mq. In other words, the 0.25 keV 
background is contributed by almost the entire range of halos cor- 
responding to groups. We assume that isolated galaxies have total 
masses of up to ~ IO^^Mq, and that an X-ray cluster with a tem- 
perature of 1 keV has a mass of ~ 10^* Mq. This result means that 
if we wish to match our upper limit on the 0.25 keV background, 
then we must reduce the contributions of almost all groups. 

The triangles in Fig. W give the mean fractional contribution 
per halo as a function of mass. The plus signs give the total number 
of hot gas halos obtained in the simulation. Therefore, the product 
of these two curves reproduce the flat plateau traced by the squares. 

The situation is quite different for the XRB at 1 keV. In Fig. ra 
the squares peak quite sharply at around 6 x 10^^ Mq. The dif- 
ference is brought about by the much steeper curve traced by the 
triangles, because less massive halos contribute less at 1 keV. 

Returning to Fig. W, we note that the sharp drop in contribu- 
tions below 4 X IO^^Mq appears to be entirely due to the low-mass 
'cutoff' in the dotted curve (cf. Fig. H). This 'cutoff results from 
the transition from cold to hot collapses, which is controlled by the 
parameter tq. Therefore, to investigate the sensitivity of the XRB 
to To, we reduced its value from 1 to 0.1. This gave the dot-dashed 
spectra in Figs. B and M which show only a small increases from 
the solid spectra (about 25 per cent at 0.25 keV). For the ACDM 
case, we show the new breakdown of contributions at 0.25 keV in 
Fig. ]n The curves now extend to a lower mass, showing that the 
plateau does indeed end around 4 x IO^^Mq , even though the tran- 




10 -^ 10 

halo total mass (Mq) 

Figure 7. As Fig. W, but for the 'W simulation in the ACDM case. The plot 
shows halo contributions to the 0.25 keV background when the luminosities 
of groups are reduced to match the extrapolated Lx~T relation from White 
et al. (1997). By doing so, the background gradually becomes dominated 
by halos at ~ IO^^Mq. In the 'PW simulation the luminosities of most 
groups are reduced much further, but the 0.25 keV background does not 
drop much as it becomes 'supported' by halos of order IO^^Mq. 



sition from cold to hot collapses has been pushed to less massive 
halos. This explains why the 0.25 keV background is not sensitive 
to the parameter tq. 

For the last two simulations listed in Table hi, we fixed the 
luminosities of halos according to observed Lx — T relations, as 
described in section [2.3| . The resulting spectra, shown in Figs. 
and W, satisfy our upper limit on the 0.25 keV background. In ad- 
dition, the 1 keV background in the OCDM case is very close 
to 1 keV cm~^ s~^ sr^^ keV^^, as expected from observations 
(see above). (In the ACDM case, the 1 keV background is slightly 
smaller.) In spite of the great difference in the Lx — T slopes 
adopted for groups, the 'W and 'PW simulations produce spec- 
tra that differ by less than a factor of 1/3. The reason for this is 
illustrated in Fig. u\ and discussed in the caption. Notice that even 
in the 'W case, the contributions of individual groups at 0.25 keV 
have been reduced by an order of magnitude or more, across the 
entire mass range. In a further simulation, we took the Lx — T re- 
lation from the 'PW simulation and 'pivoted' the slope below the 
intersection point at 0.8 keV, to give Lx oc T^ (so that luminosities 
are higher than in the 'W simulation). This increased the spectrum 
around 0.25 keV by a factor of 1/3 relative to the 'W simulation. 



3.3 Source counts 

The source counts from the above simulations are displayed in 
Figs. |g to |n} Figs. |8| and g show results from the OCDM simu- 
lations in the 0.1-0.4 keV and 0.5-2 keV bands, respectively. Like- 
wise, results for ACDM are shown in Figs, hd and hll 

As with the spectra, the log A^'-log S curves fall clearly into 
two groups in each figure, given by the first four and the last two 
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Figure 8. Simulated log A^-log S functions in the 0. 1-0.4 keV band, for the 
OCDM cosmology. The legend for the different lines is given in Table HJ 



Figure 10. Simulated log A^-log S functions in the 0.1-0.4 keV band, for 
the ACDM cosmology. The legend for the different lines is given Table In. 
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Figure 9. As Fig. W, but showing source counts in the 0.5-2 keV band 
(OCDM cosmology). The heavy dotted-line approximates the observed 
source counts. 
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Figure 11. As Fig. hol but showing source counts in the 0.5-2 keV band 
(ACDM cosmology). The heavy dotted-line approximates the observed 
source counts. 



simulations listed in Table 111 Little data exists for comparison, ex- 
cept above 10"^'* erg s~^ cm^^ in the harder band. The logA*'- 



log 5* function from the WARPS survey ( Jones et al. 1998 1 extends 



do wn to 6 X 10 ^^ er g s ^ cm ^, and the log A'^-log S function 
of Rosati et al. (1995) covers fainter fluxes down to 1 x 10"^* 



erg s cm . Both observed functions are closely approximated 
by the simple equation Af(> 5) = (10"^^erg s~^ cm"^/^), 



which lies an order of magnitude below the 'self-similar' predic- 
tions shown in Figs. S and [iTl This is to be expected, for the sources 
in this flux range are identified as small clusters in the simulation, 
with temperatures of ~ 1 keV. Hence, we already know that they 
should not be self-similar, as they lie at (or just off) the l ower end 
of the Lx - T relation of White, Jones, & Forman (1997). 



Several of the differences between the log A^'-log S functions 
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can be traced back to the spectra. For example, the dotted and solid 
curves (which differ in the metallicity used) are very close in the 
0. 1-0.4 keV band, but differ by almost a factor of 2 in the 0.5-2 
keV band. As discussed above, this can be attributed to the iron 
L complex. Reducing ro to 0.1 (dot-dashed curves) increases the 
number of low-mass halos with hot gas; this makes negligible dif- 
ference in the 0.5-2 keV band, but some change can be seen in the 
0. 1-0.4 keV band. In all cases, increasing the CDM shape param- 
eter to r = 0.25 increases the number counts by roughly 50 per 
cent, though this depends somewhat on the flux. 

The log Ai'-log S functions from the 'W and 'PW sim- 
ulations differ little from each other, especially above 10^^* 
erg s"'^ cm~^. This is because the counts are mostly dominated 
by (small) clusters in both cases. Therefore they do not depend 
on which Lx — T relation we use for groups. This suggests that 
the predictions in the 0.5-2 keV band above 10^^* erg s^^ cm"'^ 
should come close to the observed counts, since we have imposed 
the observed ix — T relation for clusters and fitted the observed 
cluster temperature function. However, the model source counts 
are a factor of 2 to 3 smaller The main reason for this discrep- 
ancy is because we do not include any scatter in the Lx — T re- 
lation used in the simulations. The XRB is not sensitive to scat- 
ter, for it measures the mean properties of halos. However, scatter 
in the flux of sources, when combined with the large number of 
faint sources compared to bright ones, can significantly increase 
the model log TV-log S function. 

We now describe a simple correction for this effect. We first 
note that the best-fit Lx — T relation of W hite, Jones, & For- 
man (1997) is calculated in logarithmic space. Therefore the sim- 
plest way to include scatter is to give logLx (at each tempera- 
ture) a gaussian distribution centred on the value given by the best- 
fit. Let the standard deviation, a, be independent of temperature 
(see their Fig. la). It follows that the new AN /AS distribution is 
given by convolving the old one with this gaussian, with log 5* as 
the independent variable. We now approximate the model log A'^- 
log S functions (for S > 10^^* erg s^^ cm^^) by the power law 
A'^(> S) oc 1/5* — as we did for the observed source counts — 
which implies AN /AS ex 1/S^. It is not hard to show that AN /AS 
remains a power law when convolved with the gaussian, but it is 
shifted upwards by a factor of 10^°^ imo jijgj-efore an intrinsic 
scatter in logLx of a = 0.3 (i.e. a factor of 2 in each direction) 
increases A'^(> S) by a factor of 2.6. This brings the model log A'^- 
log 5* functions roughly in agreement with the data. However, the 
true intrinsic scatter for small clusters is not well determined, and 
this result highlights the inherent uncertainty in modelling source 
counts. This correction should also be applied to the 0.1-0.4 keV 
band where the slope of A'^(> S) is the same. At much lower fluxes, 
a sizable contribution from groups is not ruled out (see Figs. H to 
plj ). If their scatter is very large (as discussed below), then their 
counts would be increased in the same way. 



4 FURTHER INVESTIGATIONS 

4.1 Non-gravitational heating while maintaining the same 
gas fractions 

Here we discuss an example of heating that fails remarkably to re- 
duce the simulated XRB to the required level. 

In our heating models described in WFN99, we allow the to- 
tal energy of gas halos (the sum of thermal and potential energies) 
to increase in the presence of excess energy from non-gravitational 
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Figure 12. Further .simulated spectra in the OCDM cosmology. The solid 
spectram from Fig. O is shown here for compaiison. The dashed spectram 
was obtaine d by heating all gas halos to the point of being marginally bound 
(see section p. 1[ ). For the dotted spectrum, the gas fr actio ns of halos were 
determined naturally from their progenitors (section [=t.2[ ). The dot-dashed 
spectrum was given by a preliminary study of the injection of excess energy 
by gala xies which simultaneously fits the Lx — T relations for clusters 
(section [4. 3|). 



heating. This both raises the gas temperature and flattens the den- 
sity profile of a gas halo, but we maintain the same gas mass within 
Tvir. For the isothermal gas profiles used in this paper, this is mod- 
elled by reducing the slope parameter 77. The X-ray luminosity can 
be reduced by up to an order of magnitude as a result, before the 
total energy of the gas halo becomes positive (relative to the appro- 
priate potential), at which point we regard the gas as gravitationally 
unbound. 

In two further simulations, we modify the first model in Ta- 
ble hi by heating the gas in all halos > 3 x IO^'^Mq to the point 
of being marginally bound, whilst retaining a gas fraction of 0.17. 
(100 realisations were used in these simulations.) The strong heat- 
ing roughly doubles the temperature of a halo and halves the value 
of 77. These simulations may be expected to give the maximum pos- 
sible reduction in the XRB given a gas fraction of 0.17. However, 
the resulting XRB spectra, given by the dashed curves in Figs. Il2 
and O, show a reduction of only a half (at 0.25 keV) compared to 
the solid spectrum. 

Notice that we have not heated halos below 3 x IO^^Mq . Do- 
ing this doubles the dashed spectrum at 0.25 keV, due to a whole 
new contribution from halos of ~ lO^^M© — based on the criterion 
of ro = 1, these halos become hot collapses (see section 2.1.1) 



as a result of their increased cooling times. However, we are more 
interested in how much we can reduce the XRB. 

Fig. II4 shows the new distribution of halo contributions for 
ACDM. Except for the point at 4 x lO^^M©, the total contribu- 
tions (squares) at all higher masses are reduced by the heating. (The 
contribution at 4 x IO^^Mq increases for the same reasons given 
above, but the effect on the total 0.25 keV background is small.) 
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Figure 13. As Fig. hd but for ACDM. (The solid spectrum is taken from 
Fig-i) 
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Figure 14. Plot of halo contributions to the 0.25 keV background when gas 
halos are maximally heated as described in section uA, for the ACDM case. 
The symbols are described in Fig. W 



The reason for the rather modest reduction in the 0.25 keV back- 
ground has to do with the large fraction of gas able to cool in the 
smaller groups when heating is absent. The fraction of gas that is 
able to cool decreases gradually as we progress to massive groups. 
The depletion of gas from the centres of halos can greatly reduce 
the luminosity, so that even when heating is absent the groups do 
not follow Lx oc T^ on average, but obey a steeper relation. When 



the halos are heated in the manner described, the fraction of gas 
able to cool becomes small even in small groups, and we recover 
Lx oc T^ (because t; ~ 5 in all groups) normalized to a lower lu- 
minosity. Thus, the modest reduction of the XRB may be attributed 
to the already-reduced luminosities of groups (averaged over time) 
due to cooling. This also explains why heating reduces the contri- 
butions of large groups more than for small groups, as shown by 
the sloping curve in Fig. [14 (cf. Fig. 0). 

Note that if we impose Lx oc T^ for all groups and clusters 
when heating is absent, then this corresponds to ignoring the effect 
of cooling. As we mentioned in the previous section, this results 
in a 0.25 keV background that is 5 times higher than our fiducial 
result, and therefore about 10 times higher than the dashed spectra. 

The inability of this method to reduce the 0.25 keV back- 
ground to a level that is even close to the upper limit, implies that 
the gas fractions of groups (before the gas cools) must be much 
lower than 0.17. In principle, we can remove, say, all of the gas in 
groups below 2 x IO^^Mq (see Fig. |l4), which would bring the 
0.25 keV background just below the upper limit; but the model is 
already so contrived that groups are almost certainly gas-poor up 
to ~ 1O^*M0. This conclusion is supported by results from X- 
ray studies of individual groups. Since the average gas fraction of 
clusters is at least 0.17, we argue that a large fraction of the gas 
belonging to groups must be outside their virial radii (see below). 



4.2 Gas fractions determined by inheritance 

So far, we have fixed the initial gas fractions of groups at 0.17. 
We now consider the effect of relaxing this assumption, so that the 
amount of gas in a halo is naturally determined by the amount ac- 
creted or inherited from its progenitors. We use the same param- 
eters as the first model in Table nl but modify the primordial gas 
fraction so that we obtain large clusters with gas fractions close to 
0.17. In the OCDM case, we use a primordial gas fraction of 0.25, 
and for ACDM we use 0.2 3. (F is increase d marginally by the new 
values of Q.b, according to Sugiyama 1995.) As in the fiducial sim- 
ulations, we assume no non-gravitational heating. 100 realisations 
were used in each simulation. The resulting spectra are given as 
dotted curves in Figs. Il3 and |l3[ 

The dotted spectra are surprisingly close to the fiducial results 
and almost coincide with them. This in fact hides a large scatter 
in the gas fractions of individual groups. For example, in halos of 
8 X IO^^Mq we find that initial gas fractions of 0.14 to 0.20 are 
common. Star formation was included in galaxies according the 
model described in WFN99. This consumes gas at the 10 per cent 
level (using clusters as samples of baryons), but hot gas tha t cool s 
in groups is converted into baryonic dark matter (see section [2. l.l| ). 
The effect of excess energy from supernova heating on the gas halos 
of groups and clusters is ignored in these simulations. 

Thus as far as the XRB is concerned the primary difference 
from our fiducial simulations is the freeing of the gas fractions. Al- 
though a large amount of gas is able to cool in a small group, this 
is compensated at the next collapse by new material with relatively 
high gas fractions. As a result, the average gas fraction of newly- 
collapsed groups remains close to 0.17. Intuitively, one might ex- 
pect the groups to have higher initial gas fractions than the clus- 
ters, as would be the case if clusters formed only from group-group 
mergers, but this is of course not true in reality. These results show 
that the simplification made in our fiducial simulations has remark- 
ably little effect on the predicted XRB. 

An improvement on our model would be to follow the growth 
of halos more closely in time, by refining the simulation of merger 
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trees. However, it seems that adding gas more continuously to halos 
(compared to adding it in one go when a larger block is ready to 
collapse) would be more likely to increase the predicted XRB than 
to reduce it. 



4.3 Energy injection from galaxies 



In section ll 



we found that the gas fraction in collapsing groups 
must be much lower than that in clusters, but from the discussion in 



section 12 it is clear that cooling alone cannot produce this result. 
There must be other mechanisms to prevent gas from following the 
dark matter into the halos of groups, but which do not prevent it 
from falling into rich clusters. 

Such a scenario can be achieved by giving the gas sufficient 
excess energy so that it would not be gravitationally bound to a 
group, but would succumb to the deeper potential well of a cluster. 

The minimum excess energy required to unbind a halo is 
that required to make its total energy zero relative to the appro- 
priate potential. A scatter plot of this minimum energy as a func- 
tion of halo mass is given in Fig. 7 of WFN99 (where it is called 
the 'binding energy' of the halo). However, the current work only 
demonstrates that it is necessary to drive the bulk of the gas be- 
yond ryjr and this may requir e significantly less energy (see also 
Balogh, Babul, & Patton 1999). Some clouds may be heated more 



strongly than others; while the bulk of the gas may be unbound, 
some may still collapse to form a halo with a much lower gas frac- 
tion. This complicates any estimate of the minimum excess energy 
required to explain the 0.25 keV background. However, if the bulk 
of non-gravitational heating is injected in the smaller 'branches' of 
a merger tree (;$ 5 x IO^^Mq) by the galaxies within, then the 
excess energy of gas associated with groups should be around 1~3 
keV/particle (or more to allow for dilution), since this is the excess 
energy required to fit the properties of X-ray clusters (WFN99). 
This is similar to a model proposed by Pen (1999). 

A large energy input is unlikely to be uniform, so that one 
consequence of this scenario would be a large increase in the scatter 
of the Lx — r or the Lx — <J distribution as we go below ^ 1 
keV and gas halos become (partially) unbound. An increase in the 
scatter o f these distributions is suggested by the Hickson Compact 
Groups (Ponman et al. 1996). Indeed the large scatter in properties 
also extends to the T — a distribution (see also Helsdon & Ponman 
1999). 

In the final simulations we make a preliminary study of the 
injection of excess energy in halos below 10^^ Mq. We assume that 
the energy is due to AGN and possibly supernovae. The energy is 
then retained in the gas as long as it is not radiated. In our model no 
excess energy is lost while the gas remains outside virialized halos 
(see also below). In order to simultaneously fit the constraint on the 
XRB and the Lx — T relation for clusters, we require an excess 
specific energy of around 2.8 keV/particle in cluster gas (WFN99, 
Model B for isothermal profiles only). This flattens the density pro- 
files of small clusters and reduces their luminosity. As discussed in 
our earlier paper, the excess energy in the cluster only approximates 
the actual energy injection required, which is likely to be smaller 
because a 'gravitational contribution' to the excess energy can re- 
sult when the gas is displaced by strong heating. Nevertheless, this 
is the approximation made in the simulations. 

Our prescription for energy injection is as follows: we give 
all gas that was ever associated with a halo in the range (0.015- 
l)xlO^'^M0 an excess specific energy of 7 keV/particle (in the 
OCDM case) or 5 keV/particle (ACDM). In practice, we simply set 
the excess energy of gas to this level at the end of the relevant col- 



lapse steps; this occurs even if the gas is not bound to the halo — it 
only needs to be associated with the dark matter in the halo. Due to 
dilution by unheated gas, these result in excess energies of around 3 
keV/particle in clusters, but with a significant scatter. We also lower 
the primordial gas fract ion to 0.18, but all other parameters are the 
same as in section 4.2 Since we do not model partially unbound 



halos, most groups have no gas at all in the simulations. We used 
400 realisations for each run. 

The resulting spectra are shown as dot-dashed curves in 
Figs. O and Il3. They satisfy the upper limit in both cases. More 
dilution of the excess energy seems to occur in the OCDM cos- 
mology, which therefore required a higher 'initial' excess specific 
energy. 

The resulting clusters have gas fractions of about 0.17. By 
comparing with the primordial gas fraction it is evident how lit- 
tle gas has cooled. Although we expect hardly any gas to cool in 
groups, galaxies are also strongly affected by the excess energy 
from their progenitors. In the simulations, this is exacerbated by 
the averaging of excess energies over all the gas associated with 
a collapsing halo. To avoid this problem, it may be necessary to 
inject most of the required excess energy in halos comparable to 
~ IO^^Mq. 

As a result of heating, gas would be ejected from a halo 
as a wind, terminating any star formation in the process. Such 
a wind may naturally result from the growth of a massive black 
hole (Fabian 1999). To estimate the black hole growth that is re- 
quired, let the energy available for the heating of gas be given by 
e AMbhc'^, where AMbh is the mass accreted onto the black 
hole and e is the mass-to-energy conversion rate. If we let e = 
0.1 (the value often used for the mass-to-light conversion rate), 
and distribute the energy over gas of mass Afgas, then we ob- 
tain an excess specific energy equal to O.lc^ AMbh /Afgas = 
5.8x 10''(AAfBH/Afgas) keV/particle. Thus, lO'^M© of black hole 
growth is able to supply IO^^Mq of gas with about 6 keV/particle. 
In a similar calculation in WFN99, we used the black hole den- 
sity of the universe to show that about 4 keV/particle of excess 
energy can be obtained in this way even if it is averaged over all the 
baryons in the universe. 

Once the gas is heated above a temperature of 1 keV in a 
galaxy halo, its cooling time becomes very long — comparable to 
a Hubble time, depending on its density. As the gas expands, its 
cooling time remains high (icooi oc l/(r^ " A(r, Z)) if the gas ex- 
pands adiabatically), but more importantly it converts a large frac- 
tion of its thermal energy into potential and kinetic energies. In this 
way, the gas should be able to retain most of its excess energy until 
it recoUapses into a cluster. 



5 IMPLICATIONS FOR THE IGM AND GALAXY 
FORMATION 

In the previous section, we first showed that the group population as 
a whole has a much lower gas fraction than X-ray clusters, and then 
proposed non-gravitational heating at the level of ~ 1 keV/particle 
as a means of accounting for both the properties of groups and X- 
ray clusters, in a natural and self-consistent manner. Indeed, the low 
gas fractions of groups and the large intrinsic scatter in their proper- 
ties compared to clusters (section 4.3) both lend support to this high 



level of heating, whi ch we originally proposed for clusters. Simi- 
larly, Renzini (1997) noted a precipitous drop in the intra-cluster 
medium (ICM) mass-to-light ratio at around 1 keV, as well as in 
the ICM iron mass-to-light ratio at the same temperature ('light' 
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refers to the total B-band luminosity of the galaxies in the cluster). 
He concluded that clusters could not have formed by assembling 
groups similar to those observed. In the heating scenario, this is 
resolved by allowing the gas to 'reunite' with the groups on the 
formation of a cluster. 

In the proposed model, the high energy of heated gas in the 
IGM means that it cannot remain in the filamentary and sheet-like 
structures seen in N-body simulations. The actual distribution of 
the different gas phases in the IGM would depend on details of how 
heating occurred, a problem analogous to supernova heating in the 
inter-stellar r nedium. This is contra ry to the evolution of baryons 
described by ^en & Ostriker (1999). In their cosmological hydro- 
dynamic simulations, almost half of the baryons at redshift zero 
lie in the temperature range 0.01-1 keV and exist in filaments and 
more clumpy structures. This gas is heated primarily by collisions, 
where the energy has a gravitational origin. When heating is in- 
cluded at the proposed level, a large fraction of this gas would be 
expelled from the potential wells of filaments and groups, result- 
ing in a smoother and more diffuse gas distribution. As a result, the 
likelihood of detecting this la rge reservoir of baryons in the uni - 
verse would be much reduced (Fukugita, Hogan, & Peebles 1998). 



Since even the smallest groups are affected by the heating, it 
seems likely that the heating process would play an important part 
in the evolution of galaxies. For example, the high end of the lumi- 
nosity function of galaxies (where the Schechter function decays 
exponentially) may be influenced by the heating process. Assum- 
ing that the source of energy is AGN, the possibility that AGN and 
galaxy formation are intimately connected will continue to be a 
topic of much discussion. 



6 SUMMARY OF CONCLUSIONS 

We have made a systematic study of the soft X-ray background 
from hot gas halos, with the aim of constraining the properties of 
the gas halos of groups. Using Monte Carlo simulations of halo 
merger trees coupled with realistic gas density profiles, and includ- 
ing the effects of gas removal due to cooling, we calculated the 
XRB spectrum along with source counts in the 0.1-0.4 keV and 
0.5-2 keV bands. In addition, we investigated the composition of 
the XRB in terms of the masses of groups that contribute at 0.25 
keV and 1 keV. Our main conclusions are as follows: 

• Radio-quiet quasars are able to account for almost all of the 
extragalactic XRB at 0.25 keV. As a result, we set an upper limit of 
4 keV cm~^ s~^ sr~^ keV~^ on the 0.25 keV background due to 
gas halos. 

• In the absence of non-gravitational heating, the predicted 0.25 
keV background is an order of magnitude higher than this upper 
limit. In addition, it is contributed by the entire mass range of 
groups, from ~ 5 x 10^^ to ~ IO^Mq. 

• The removal of gas due to cooling plays an important part in 
determining the XRB in this case. Excluding this effect increases 
the predicted 0.25 keV background by about a factor of 5. 

• Maximally heating the gas halos of groups without changing 
their gas fraction reduces the 0.25 keV background by only a factor 
of two or less. It follows that most of the gas associated with groups, 
down to halos of ~ 5 x lO^^M©, must be outside the virial radii of 
these halos. 

• The properties of both groups and X-ray clusters can be natu- 
rally explained by a model in which the gas is given (on average) 
excess specific energies of 1-3 keV/particle by non-gravitational 
heating. 



• In addition to satisfying the constraint on the XRB, this would 
result in a large scatter in the X-ray properties of groups (assuming 
that the heating is inhomogeneous), as well as a dichotomy in the 
properties of gas halos above or below T ~ 1 keV, both of which 
are supported by observations. 

• The excess energy is most likely injected by galaxies in the 
smaller 'branches' of the halo merger tree (;$ 5 x IO^^Mq), by 
active galactic nuclei and possibly supernovae. 

• This greatly reduces the likelihood of detecting the large reser- 
voir of cosmic baryons which would otherwise be expected in 
groups and filaments. 
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